Stochastic force generation by small ensembles of myosin II motors 
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Forces in the actin cytoskeleton are generated by small groups of non-processive myosin II motors 
for which stochastic effects are highly relevant. Using a crossbridge model with the assumptions 
of fast powerstroke kinetics and equal load sharing between equivalent states, we derive a one-step 
master equation for the activity of a finite-sized ensemble of mechanically coupled myosin II motors. 
For constant external load, this approach yields analytical results for duty ratio and force-velocity 
relation as a function of ensemble size. We find that stochastic effects cannot be neglected for ensem- 
ble sizes below 15. The one-step master equation can be used also for efficient computer simulations 
with linear elastic external load and reveals the sequence of build-up of force and ensemble rupture 
that is characteristic for reconstituted actomyosin contractility. 

PACS numbers: 87.10. Mn.87.16.Ln,87.16.Nn,82.39.-k 



Generation of motion and force by ATP-powered 
molecular motors is a hallmark of living systems [l[ . In 
their cellular environment, molecular motors usually op- 
erate in groups 0. A striking example is force genera- 
tion in skeletal muscle, where hundreds of non-processive 
myosin II motors are assembled into the thick filaments of 
the sarcomeres. Since the pioneering work of Huxley Q, 
the statistical physics of large ensembles of myosin II mo- 
tors has been studied in great detail. It has been shown 
that in order to describe the response of skeletal muscle 
to varying loading conditions, it is essential that the un- 
binding rate of myosin II from actin is strain-dependent 
and decreases under load 0, 0]- In contrast to e.g. the 
processive motor kinesin, this makes myosin II a catch 
rather than a slip bond p, 0] and leads to recruitment of 
additional crossbridges under load Q. 

The collective activity of myosin II motor ensembles is 
also essential for the generation of motion and force in the 
actin cytoskeleton of non-muscle cells. In this case, the 
actin structures are far more disordered than in muscle 
and non- muscle myosin II is usually organized in mini- 
filaments comprising 10 — 30 motors Q. For such small 
numbers of motors, stochastic effects will become impor- 
tant and are indeed observed in experiments. Measure- 
ments of tension generated by myosin II motors in recon- 
stituted assays, e.g. in three bead assays [6, 10, llj, active 
gels [H, or motility assays Q, reveal noisy trajecto- 
ries, typically with a gradual increase of tension followed 
by an abrupt release, which is likely due to detachment 
of the whole ensemble {slip). However, a detailed and 
analytically tractable description for this biologically im- 
portant situation is still missing. 

The collective activity of mechanically coupled molecu- 
lar motors has been investigated before in the framework 
of a generic two-state Fokker-Planck equation in which 
ensemble size enters into the noise intensity 14, 1^. In 
order to study effects of molecular details for ensembles 
of myosin II motors, crossbridge models originally devel- 



oped for skeletal muscle can be used as a starting point 
[j, 0] . Due to their complexity, these models are usually 
studied by computer simulations. Analytical progress has 
been made with a mean field approximation for large sys- 
tem size [l6j . Exploiting a separation of time scales in the 
myosin II cycle and using the assumption of equal load 
sharing between motors in equivalent states, here we de- 
rive a one-step master equation which explicitly includes 
the effects of catch bonding and small system size. A 
one-step master equation has been introduced before for 
transport by finite-sized ensembles of processive motors 
with slip bond behavior , but not for non-processive 
motors with catch bond behavior. Our results suggest 
that stochastic effects are particularly important for en- 
semble sizes below 15, which corresponds to the typical 
size of cytoskeletal mini-filaments. 

Model. We model the myosin II cycle by three discrete 
mechano-chemical states. The cycle is shown schemat- 
ically in Fig. [1^. To allow for comparison with earlier 
work, transition rates and most other molecular parame- 
ters are taken from Refs. 0, Hj. In practice, they will 
depend on ATP concentration and the exact type of 
myosin II @, [l3|- In the unbound state (0), the motor- 
head is loaded with ADP and Pi and the lever-arm is 
in its primed conformation. The motor then reversibly 
transitions to the weakly bound state (1) with forward 
rate /cqi — 40 s~^ and reverse rate fcio — 2s^^. After re- 
lease of Pi; the lever-arm swings to the stretched confor- 
mation and the motor enters the post-powerstroke state 
(2). The transition rates between the two bound states 



are relatively high, with fci2 ~ fc2i 



10^! 



Replacing 



ADP by ATP, unbinding from the substrate and hydrol- 
ysis of ATP brings the motor back to the unbound state 
(0). This last step is irreversible, with rate k2o — 80 s~^. 
Most important in our context, both powerstroke and 
unbinding depend on load. The powerstroke (1) — >■ (2) 
moves the lever-arm forward by d ~ 8 nm and strains 
the elastic neck-linker. Unbinding from (2) requires fur- 
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FIG. 1. (a) Myosin II motor cycle with three mechano- 
chemical states, (b) Mechanical coupling of a motor ensemble 
moving to the right. The external load pulling to the left is 
balanced by elastic forces in the neck-linkers of the motors. 
Swinging of the lever-arm increases the strain of the neck- 
linker and hence the force exerted by a motor, (c) In the par- 
allel cluster model (PCM), all motors in the same mechano- 
chemical state have the same strain, (d) The correspond- 
ing two-dimensional reaction network, with irreversible tran- 
sitions marked by arrows, (e) The effective one-dimensional 
network following from the assumption of local thermal equi- 
librium (LTE) of the bound states. 



ther movement of the lever-arm, thus making unbinding 
slower under load (catch bonding). 

As shown schematically in Fig. [T|d, the (upper) motor- 
filament mechanically couples the different motors in an 
ensemble to each other. Due to the strain-dependence 
of the rates, they are also dynamically coupled. Hence, 
a complete description of the ensemble dynamics has to 
include conformational state and strain of every motor. 
To arrive at a tractable model, we first note that the 
motors pull in parallel. We next assume that all motors 
in the same mechano-chemical state exert the same force 
and hence have the same neck-linker strain. We thus 
arrive at the parallel cluster model (PCM) depicted in 
Fig. [TJ;, in which the state of an ensemble with A^t motors 
is characterized by the number i of bound motors and the 
number j < i of motors in the post-powerstroke state. 
The number of motors in the weakly bound state follows 
as i — j. In the PCM, each motor in the weakly bound 
state has the same strain Xij, where the indices indicate 
the dependence of the motor strain on the ensemble state 
(i, j). The powerstroke stretches the elastic neck- linker 
by d, so that motors in the post-powerstroke state have 
the strain Xij -\- d. The strain Xij of the weakly bound 



motors follows from the balance of the external load Fcxt 
and the elastic motor forces: F^xt = kin[ii—j)xij-\-j{xij-\- 
d)]. Here km — 2.5pNnm~^ is the spring constant of the 
neck-linkers. For Fcxt = const, the force balance leads to 

Xij = (Foxt - jkmd) /ikm ■ (1) 

Thus the strain Xij of the weakly bound motors is a 
state variable determined by external load and both bind- 
ing and powerstroke dynamics. If all motors are in the 
weakly bond state (j = 0), it is positive. It can become 
negative if sufficiently many motors have gone through 
the powerstroke and if the external load is not too large. 
The strain Xij -|- d of the post-powerstroke motors always 
stays positive and eventually drives force generation and 
motion. 

In the PCM, the network of reactions between states 
is two-dimensional (see Fig. [TJi). Due to slow 
binding and unbinding, local thermal equilibrium (LTE) 
is maintained for the bound states [5|. When i mo- 
tors are bound, the probability that j motors are in 
the post-powerstroke state follows the Boltzmann dis- 
tribution p{j\i) = exp{—Eij/kBT)/Z, where Z is the 
appropriate partition sum. The energy Eij = E^i -\- 
jEpp + Ecxt in state (i, j) is the sum of elastic energy 
Eel = fcm((i — + j{xij -\- d)^)/2 stored in the neck- 
linkers, free energy bias Epp ~ — 60pNnm towards the 
post-powerstroke state, and a possible external energy 
contribution E^xt- For F^xt = const, we have -Ecxt = 0. 

LTE of the bound states allows us to project the j- 
axis onto the i-axis, thus arriving at a one-dimensional 
reaction scheme with index i as shown in Fig. [1^. Then 
the probability Pi{t) that i motors are bound at time t 
obeys the one-step master equation 

^p^ = r{i + l)pi+i + g{i - l)p^-l - [r{i) + g{i)]pt . (2) 

The probability to find an ensemble in state is 
Pij(t) = p{.j\^)Pi{t) ■ The forward rate is g(z) = {Nt—i)kQi 
because Nt — i free motors can bind. Unbinding is pos- 
sible from states (1) and (2) so that the reverse rate for 
given i and j is r(i,j) = {i - j)fcio + jk2o{i,j)- Aver- 
aging over j gives r(i) ~ pij\i)'>'ii, j)- The off-rate 
from state (2) depends on the applied load as k2oi'i,j) = 
ko cx.p {-km{xij -\- d)/Fo), where Fq ~ 12.6pN. The 
strain-dependence of fc2o makes myosin II a catch bond. 
With these prescriptions, the one-step master equation 
Eq. ^ is fully specified for the case of constant external 
load, Fcxt = const. If the external load depends on the 
position of the ensemble, like in the case of linear elastic 
loading, Eq. ^ has to be solved together with additional 
prescriptions for ensemble movement (sec below). 

Binding dynamics for constant load. Mathematically, 
the reduction to Eq. ([2]) is a dramatic simplification, be- 
cause many general results are known for one-step master 
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FIG. 2. (Color online) Analytical results for the parallel clus- 
ter model (lines) and computer simulations with individual 
motor strains (symbols) for constant external load, (a) Av- 
erage number A^b of bound motors following from Eq. Q as 
function of Fext/At for At = 4, 8 and 15. (b) Duty ratio pd 
given by Eq. © as function of At for Fext/At = 0.013 pN, 
1.262 pN, 3.786 pN and 8.834 pN. 



equations jl8j . The stationary distribution is 

llj=0 r( 



(3) 



i=0 r-O-l-l) 



Fig. Hi plots the average number of bound motors, A'b = 
if) ~ J2i=o'''Pii'^)i function of Fext for different en- 
semble sizes Nt (lines). The increase of A^b is due to the 
catch bond character of the post-powerstroke state. With 
increasing load, r{i,j) decreases to {i — j)kio, so that r(i) 
is small because p{j\i) is strongly biased towards large j 
and fcio ^ ^20- For skeletal muscle, the recruitment of 
additional crossbridges under load has been observed ex- 
perimentally as predicted by computer simulations Q . 
Here it follows in a relatively simple way from analyti- 
cal considerations. In order to validate the PCM leading 
to Eq. ^ , in Fig. [5^ we also show results of computer 
simulations which incorporate an individual strain value 
for each motor (symbols). The agreement is very good, 
except at very small load, where differences in the strain 
values between different motors reduces unbinding, an 
effect which is less relevant under larger load. 

Next we discuss the effect of system size Nt. In general, 
smaller ensembles are more likely to detach as a whole. 
The mean first passage time for ensemble detachment 
after binding of the first motor is 



Nt 



h '^^^ ii 



(4) 



It is a polynomial of order At — 1 in the ratio of binding to 
unbinding rate and increases exponentially with ensemble 
size. Once the ensemble has detached, on average it takes 
the time Tqi = l/g{0) = 1/Ntkoi to rcbind. Wc define 
the duty ratio of an ensemble as 



Pd = Tio/ (Tio + Tqi) 



(5) 



Fig.Hb plots Pd as function of A't for different Fext (hues). 
Because Tio increases and Tqi decreases with A't, the 



duty ratio increases quickly with A't and reaches unity for 
ensemble sizes around A^t — 15- With increasing force, pd 
increases faster because of the increasing N\,. Again the 
agreement with the simulation of the crossbridge model 
with individual motor strains (symbols) is rather good 
except at very small force. Stochastic effects are expected 
to be important for duty ratios below unity, i.e. below 
ensemble sizes around 15. This implies that myosin mini- 
filaments in the cytoskeleton arc typically at the verge of 
stochastic instability. 

Ensemble movement. We now consider the spatial co- 
ordination schematically depicted in Fig. [T]d, that is, we 
assume an immobile substrate over which an ensemble 
moves to the right. The PCM assumes that all bound 
motor-heads are at the same position, which we denote 
by the coordinate z. The anchors of the motors in the 
(upper) motor-filament are located at the common posi- 
tion z — Xij . Note that whereas z increases to the right in 
Fig. [T)d, external load Foxt and strain of the motors are 
defined in the opposite direction. For the case of con- 
stant external load, F^xt = const, z is a variable which 
is slaved to the binding dynamics. When the ensemble 
works against a linear external load, i^ext = kf{z — Xij), 
the value of z enters the force balance and hence feeds 
back into the system state. In addition, here one has to 
include an external elastic energy E^xt — k[{z — Xy )^/2, 
where kf is the external spring constant. 

Although more complicated assumptions might be pos- 
sible, here we make the following simple assumptions 
for the dynamics of z within the PCM. Starting with 
a state we assume that a motor binds to the sub- 

strate with vanishing strain at the position z — Xij just 
below its anchor point. Binding of a new motor thus 
changes the average position z of the bound motor-heads 
by Az = {iz + {z — Xij)) /{i + I) — z = -Xij/{i + 1). 
Unbinding docs not change z, because all bound motor- 
heads are at the same position. The powerstroke does 
not change z either, because it does not affect the posi- 
tions of the motor-heads. However, Xij is affected by 
binding and unbinding as well as by the powerstroke 
via the force balance, so that the position of the motor- 



filament z - 



is affected by all these transitions. When 



an ensemble detaches completely from the substrate, the 
motor-heads relax to the position z — Xij of the anchors. 
The detached ensemble then moves backwards with ve- 
locity Vs = — TyJ^ext (slip), where 77 is the effective mobility 
of the motor-filament. 

With these additional prescriptions, the rates defined 
for Eq. ([2]) can now be used to investigate the details of 
the stochastic movement of the motor ensemble for arbi- 
trary laws for the external load. To simulate stochastic 
trajectories, we use the Gillespie algorithm [l^. After 
every change of i, Xij and p{j\i) are updated to cal- 
culate the average strain of the weakly bound motors, 
Xi = jyj^QXijp{j\i), and the transition rates r(i) and 
g{i). In case of binding, we change z by Az = —Xi/{i+l). 
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FIG. 3. Stochastic trajectories, (a) Constant load: Average 
head position z (upper panel) and number i of bound motors 
(lower panel) as function of time t for ensemble size A^t = 8 
and load Kixt = 1.26pNA't. (b) Linear load: Average head 
position z (upper panel) and number i of bound motors (lower 
panel) as function of t (Nt = 4, kt/Nt = 0.0126 pN nm~^). 
In (a) and (b), an detached ensemble slides backwards with 
mobility rj = 10^ nmpN"^ s^^. 



Fig. 13^ shows a stochastic trajectory of an ensemble 
working against constant load. The lower panel shows 
the number of bound motors i, the upper panel the av- 
erage head position z as function of time. When bound, 
the ensemble moves forward with fluctuations around a 
steady state velocity. A slip leads to backsteps of aver- 
age size VsTqi. Fig. [3}3 shows a trajectory for an ensem- 
ble working against a linear elastic load. The ensemble 
is slowed down by the load building up by the forward 
motion. An increasing load stabilizes the ensemble be- 
cause A^b increases. However, the very small ensemble 
frequently detaches before reaching the stall force. De- 
tachment leads to a noisy trajectory in which the load 
fluctuates around an effective stall force. This type of 
trajectories, with gradual buildup and quick release of 
tension, resembles those experimentally observed in three 
bead assa 
assays 



ssays 

s 



s [a, |10|, lUj , active gels [12|, |13| and motility 



Force-velocity relation for constant load. In state 
(i,j), one can identify the ensemble velocity with Vij = 
—g(i)xij / {i + 1) (the ensemble only moves to the right 
when the strain defined to the left is negative). The av- 
erage stationary velocity of a bound ensemble is 



Nt 



^^v,jp{j\i)p,{oo) 



(6) 



i=l 3=0 



with Pi(oo) from Eq. ([3]). This is the force- velocity rela- 
tion of the bound ensemble at constant load. Fig.|3^ plots 
Wb as function of the external load per motor for differ- 
ent Nt. With increasing load, the velocity decreases. The 
upward convex shape of Wb(^ext) is due to the increase of 
A^b with Fext, which allows the ensemble to resist larger 
forces. For small ensembles for Fe^t/Nt > 0, bound ve- 
locity Vh and also the stall force increase with increasing 
A't. Above A^t — 15, the force- velocity curve is indepen- 
dent of A^t . This confirms our conclusion from the duty 
ratio that stochastic effects cannot be neglected up to a 



FIG. 4. (Color online) Force-velocity relation for constant 
load. Main panels: Analytical results for (a) the average 
bound velocity Vh (see Eq. ((B))) and (b) the average effective 
velocity Vea (see Eq. ((Tjl) of an ensemble as function of the 
external load per motor Fcxt/Nt for At = 4, 8, 15, 25 and 
50. For Ucff the free mobility is 77 = 10^ nmpN~^ s~^. Insets: 
Comparison of analytical results (lines) with computer simu- 
lations of the crossbridge model with individual motor strains 
(symbols) for (a) Wb and (b) Wcfi as function of external load 
Fcxt for the same parameters as in the main panels. 



system size of 15 (compare Fig. [^t). 

Assuming that the stationary velocity is established 
quickly after binding to a substrate, the walk-length of 
a motor-filament in one attachment event is given by 
dw = I'bTio- Although the bound velocity Wb decreases, 
increases with i^ext because the detachment time Tio 
increases strongly. Only upon passing the stall force, the 
walk-length drops to negative values. Comparison with 
numerical solutions of the master equation (not shown) 
reveals that dw = VhTio is a good approximation except 
for very small values of c?w Because the sliding velocity 
is negative, Vs < 0, the effective velocity is reduced by 
the occurrence of slip events: 



Vcff 



VhTio + VsToi 



Toi 



< Vh ■ 



(7) 



Fig. Hb plots the effective velocity Vf.ff as function of the 
external load per motor. Because the duty-ratio increases 
and the rebinding time decreases with A't, i.e., detach- 
ment is less frequent and backsteps are smaller, the ve- 
locity at small Fext now increases with A't. In addition, 
detachment of small ensembles leads to a faster decrease 
of Weft under load and a smaller stall force. Moreover, 
detachment leads to large fluctuations of z at the effec- 
tive stall force: instead of being stationary, the ensemble 
alternates between slow forward motion when bound and 
fast backward slipping when detached. Above the thresh- 
old of A't — 15, where the duty ratio is close to unity, the 
effective velocity is identical to the bound velocity. 

The insets in Fig. 2] compare the analytical results us- 
ing the PCM to the computer simulations without PCM 
for Vh and Veff as function of i^oxt- The agreement is 
rather good. Due to the molecular friction resulting from 
differences in strain, the bound velocity at vanishing load 
now decreases with A't and the curvature of the force- 
velocity relation is less pronounced. The stall force and 
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the role of ensemble size for stochastic effects is predicted 
well. 



Discussion. In this Letter, we have derived a mathe- 
matically tractable model for the collective behavior of 
small ensembles of myosin II motors as a function of sys- 
tem size. Our main assumption, the parallel cluster model 
(PCM) for the load sharing, was validated by computer 
simulations of a cross-bridge model with individual motor 
strains. These assumptions decrease the disorder in the 
motor strains, so that the model cannot describe power- 
stroke synchronization through load as it has been done 
before with a detailed model for skeletal muscle Q . How- 
ever, our model makes accurate predictions for central 
quantities such as duty ratio and force-velocity relation 
as a function of ensemble size. For processive motors, 
the strains of the motors are homogenized because fast 
moving motors are slowed down by the increasing load. 
For the non-processive motors studied here, this mecha- 
nism cannot operate. However, here the differences in the 
strain of the bound motors are reduced by the small duty 
ratio, thereby making the PCM a reasonable assumption 
for our purposes. Due to its computational simplicity, 
in the future the approach introduced here can be used 
for studies of the intriguing interplay between actin fil- 
aments and small ensembles of myosin II motors in the 
actin cytoskeleton of non-muscle cells and reconstituted 
actomyosin systems. 
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